close all
clear all
clc

delta=0.083;
tauI=0.0549;

i=1;

for taud=0:0.01:1
    
    Theta(i)=delta/((1-taud)*(1+tauI));
    
    TauD(i)=taud;
    
      i=i+1;
      
end


figure(1)
plot(Theta,TauD)
xlabel('Borrowing Limit')
ylabel('Dividend Tax')
xlim([0 0.5])

figure(2)
AA=area(Theta,TauD)
AA(1).FaceColor = [0.81 0.81 0.81];
xlabel('Borrowing Limit (theta)')
ylabel('Dividend Tax')
xlim([0 0.5])
